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We present temperature dependent magnetic neutron diffraction measurements on 
Ba(Fei_ a: Co a ;)2As2 for x = 0.039, 0.022, and 0.021 as-grown single crystals. Our investiga- 
tions probe the behavior near the magnetic tricritical point in the {x,T) plane, Xt r ~ 0.022, as 
well as systematically exploring the character of the magnetic phase transition across a range of 
doping values. All samples show long range antiferromagnetic order that may be described near 
the transition by simple power laws, with (3 = 0.306±0.060 for x = 0.039, /3 = 0.208±0.005 for 
x = 0.022, and = 0.198±0.009 for x = 0.021. For the x = 0.039 sample, the data are reasonably 
well described by the order parameter exponent j3 = 0.326 expected for a 3D Ising model while 
the x — 0.022 and x — 0.021 samples are near the /3 = 0.25 value for a tricritical system in the 
mean-field approximation. These results are discussed in the context of existing experimental work 
and theoretical predictions. 

PACS numbers: 74.70.Xa, 74.62. Dh, 75.50.Ee, 75.40.Cx 



I. INTRODUCTION 

Superconductivity has been a major research interest 
of the scientific community ever since the first set of ex- 
periments in 1911 suggested electron conductance with- 
out resistance, and the classes of materials that show a 
superconducting state have grown extensively over the 
yearsP From a technological standpoint, high tempera- 
ture (high- To) superconductors are particularly attrac- 
tive, such that when record breaking superconductivity 
was reported in 1986 for copper oxide based materials^ 
the field experienced an enormous surge in activityP Re- 
cently, in 2008,^ a new paradigm was discovered when 
iron pnictides were shown to display superconductivity 
at temperatures in the 50 K ranged Both the copper ox- 
ide and the iron pnictide systems are characterized by 
competition between antiferromagnetism and supercon- 
ductivity. 

With an eye to understanding better the underlying 
physics in high-Tc materials, the iron pnictides quickly 
became the subject of intense scrutiny. As a result, a 
wide variety of iron pnictides, oxypnictides, and chalco- 
genides have been found to exhibit superconductivity, 
and it has become common parlance to refer to differ- 
ent structural classes by listing the subscripts of their 
undoped chemical formulae as in '1111' for compounds 
isostructural to LaFeAsO, '122' for compounds isostruc- 
tural to BaFe2As2, or even more complicated formulae. 5 
It is starting to become clear which classes of iron based 
compounds show technological promise; it is worth men- 
tioning that while cuprates presently retain superiority in 
terms of high transition temperatures and critical fields, 
the mechanical properties of some iron based materi- 
als which are metallic may better lend themselves to 
wire production than the more brittle, layered ceramic 
cuprates. 



The interplay of superconductivity and magnetism 
seen in cuprates is also seen in the iron pnictides, which 
has motivated concentrated study of the role of the mag- 
netism in this new class of high temperature supercon- 
ductors. Interestingly, there is compelling evidence that 
superconductivity can coexist with long range iron anti- 
ferromagnetic order in the iron superconductors, while 
it remains uncertain if copper magnetic order (other 
than short range) coexists with superconductivity in the 
cupratesPThis coexistence has been shown in some '122' 
and '1111' iron compounds, where long-range magnetic 
order and superconductivity arise from bands derived 
from iron 3d electrons, but is not generic to the system 
as some formulations do not show coexistenceP For the 
examples without a sharp boundary between supercon- 
ducting and magnetic phases, the interaction and compe- 
tition between the phases is seen as a reduction in the or- 
dered moment at the onset of superconductivity.^ Other 
non-high- Tq magnetic superconductors have shown co- 
existence, such as borocarbides (i?Ni 2 B 2 CJPl, the Chevrel 
phases (iJMoeSg^, and ruthenates (RuSr 2 GdCu 2 8 ;P, 
but the most striking analogy is drawn when the same 
electronic bands par ticipa te in both phases such as in 
UPt^ and UNi 2 Al^El. Additionally, providing even 
more evidence for the importance of magnetism in the 
high-T<7 iron based superconductors is the observation of 
a magnetic resonance in the inelastic neutron scattering 
spectrum, just as previously seen in cuprates and heavy 
fermion materials 

In the present study, we investigate the cob alt-do ped 
'122' system, Ba(Fei_ x Co 2; ) 2 As 2 , at low dopingP^ESl The 
phase diagram, reproduced using measurements from 
reference^ in Fig. l] shows a typical response to dop- 
ing. The parent BaFe 2 As 2 phase shows two phase transi- 
tions when cooling below room temperature, a structural 
transformation from tetragonal to orthorhombic symme- 
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FIG. 1. (color online) Ba(Fei_^Co a; )2 As2 phase diagram in 
the vicinity of the tricritical point. The three phases in this 
region are delineated: the paramagnetic (PM) to antiferro- 
magnetic (AFM) transition at the Nel temperature (TW), the 
tetragonal (I4/mmm) to orthorhombic (Fmmm) transition 
at the structural transition temperature (Ts), and the metal- 
lic (M) to superconducting (SC) transition at the critical tem- 
perature (Tc)- For the magnetic transition, open squares 
correspond to a first order transition and closed squares con- 
nected by a solid line correspond to a second order transition. 
The vertical dashed lines show the doping levels studied in 
this work. The data points for the phase lines are from mag- 
netic susceptibility studies reported in reference^. 



try (whose temperature is denoted Ts), and a magnetic 
transformation from paramagnetic to antiferromagnetic 
(whose temperature is denoted Tjv)P^ While nearly con- 
current in the parent, doping causes a progressive separa- 
tion of Ts and Tjy with increasing x. Furthermore, upon 
initial doping, the structural transition is 2nd order and 
the magnetic transition is 1st order, and, at a doping 
denoted xtr, near the region of onset of superconductiv- 
ity, the structural transition remains 2nd order while the 
magnetic transition crosses over from 1st to 2nd orderP^S 
This crossover point in the (x,T) plane, which is near 
T = 100 K and x = 0.022, is almost certainly a tricritical 
point P 1 

To explore the nature of the magnetic phase transi- 
tions, we have performed temperature dependent neu- 
tron diffraction measurements for samples whose Co 
concentrations, x, are in the vicinity of the tricritical 
value, x = 0.022. Due to the antiferromagnetic na- 
ture of iron based superconductors, neutron scattering 
has proven ideal for studies of the magnetic structures 
and excitations In Sec. II, we outline our experimen- 
tal procedures for sample preparation and spectrometer 
configuration. Section III shows our diffraction data with 
model fits, while Sec. IV discusses the results of the fits 
in detail. Finally, in Sec. V we give our final conclusions 
and summarize the results. 



II. EXPERIMENTAL PROCEDURE 

A. Synthesis 

Single crystals of Ba(Fei_ a; Co ;c )2As2 with cobalt dop- 
ing values, x, of 0.039, 0.022, and 0.021 were grown us- 
ing a self-flux method, with details available in a pre- 
vious report!^ Samples x = 0.021 of mass 117.2 mg 
(T N w 105 K) and x = 0.039 of mass 81.7 mg (T N « 66 K; 
Tc ~ 11 K) are the samples previously used in mapping 
the phase diagram for cobalt doping 17 , and x — 0.022 of 
mass 114.6 mg (TV ~ 100 K) is from the same batch of 
the x — 0.022 material in the same study. 

B. Instrumentation 

Neutron diffraction experiments were performed on the 
BT-7 thermal triple-axis spectrometer at the NIST Cen- 
ter for Neutron Research}^ with a collimation of open- 
50'-sample-50'-120'. The x = 0.022 sample was also mea- 
sured on the high-resolution SPINS cold-source triple- 
axis spectrometer with a collimation of open-80'-samplc- 
80'-open. Both machines use the (0 2) reflection of 
pyrolytic graphite (PC) as a monochromator and ana- 
lyzer. PC filters to reduce higher order neutrons were 
employed on the BT-7 spectrometer using a fixed neu- 
tron energy of 14.7 meV (A = 2.36 A), and SPINS used 
a fixed neutron energy of 5.0 meV (A = 4.05 A) with a 
cold Be filter. Samples were mounted in the (H L)o 
scattering plane and placed inside a helium flow cryo- 
stat, and temperature control was performed in a cali- 
brated geometry capable of at least 50 mK stability. The 
energy resolution on BT-7 in this configuration is ap- 
proximately 1 meV and the energy resolution on SPINS 
in this configuration is approximately 0.2 meV. Resolu- 
tion corrections to the intensity were performed using 
the Cooper-Nathans approximationP3 Tabulated values 
for scattering lengths^ and magnetic form factors were 
usedPl 



III. NEUTRON DIFFRACTION 

To begin, we cooled each sample to less than 10 K 
and performed 8-29 scans of the (004) o , (200)o, 
(202)o, and (204)o nuclear reflections as well as the 
(103)o magnetic reflection, where the subscript denotes 
orthorhombic notation (typical lattice parameters of 
a w 5.62 A, b w 5.57 A, c « 12.94 A ). Typical data are 
shown in Fig. [2] The scale-factor was determined from 
the nuclear peaks using the BaFe2As2 Fmmm (space 
group No. 69) structure,^ with the appropriate substi- 
tutions of cobalt for iron in the structure factor cal- 
culation. In these doping ranges, it has been shown 
that the magnetic structure remains commensurate as 
in the parent phased In this way, we were able to ex- 
tract the size of the low-temperature ordered moment 
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to be 0.49 ^ B ±0.01 [i B (x = 0.021), 0.25 ^ B ±0.01 [i B 
(x = 0.022), and 0.31 /i B ±0.02 fi B (x = 0.039). 

In order to understand the critical behavior of the sub- 
lattice magnetization as a function of cobalt doping, the 
temperature dependence of the magnetic (103)o peak 
was measured on warming through the magnetic tran- 
sition for all three samples using 6 — 28 scans and inte- 
grating the intensity. Near the transition temperature, 
the intensity, I, of the magnetic diffraction peak may be 
fit to a simple power law, 
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where Tjy is the Nel temperature, A is a proportionality 
constant, and (3 is the order parameter critical exponent. 
For doped samples like those that we are studying, in- 
homogeneities in the growth process may give rise to a 
distribution of Nel temperatures, which if assumed to be 
Gaussian adds an additional term to eq. 1 for the stan- 
dard deviation, a, such that, for T < Tjy, 
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Finally, renormalization group techniques show that 
logarithmic corrections lower the effective exponents that 
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FIG. 2. (color online) 8 — 26 scans in the {H L)o scattering 
plane for Ba(Fei_ a: Co^ As2. The (103)o magnetic reflection 
present at the lowest angle is shown amplified by 50 times 
for clarity. The slightly higher background for the x — 0.039 
sample is due to a larger detector arm distance used for that 
measurement. Uncertainty bars are smaller than the data 
points, and the lines are results of model fits that are used to 
extract overall scale factors. 



one measures when approaching a tricritical point such 
that the intensity might be modeled as 
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or for the case of distribution of Neel temperatures, for 
T < T N , 
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Results of the measured temperature dependence along 
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FIG. 3. (Color online) Intensity of the magnetic (103)o peak 
upon warming. Experimental data are represented by trian- 
gles for x = 0.021, circles for x — 0.022 with both thermal 
(thick black circles) and cold (thin red circles) neutron diffrac- 
tion shown, and squares for x = 0.039. The results of model 
fits are the lines overlaying the data points, details of which 
are described in the text and in Table I. The uncertainty bars 
are derived from counting statistics and represent one stan- 
dard deviation. 
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TABLE I. Results of fitting to Eq.'s 1-4. 
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FIG. 4. (Color online) Intensity of the magnetic (103)o 
peak upon warming on a log-log scale. The three samples, 
x = 0.021 (green triangles), x — 0.022 (red circles), and 
x — 0.039 (black squares) are shown here normalized to their 
fit functions to illustrate the differences in /beta. The results 
of model fits to Eq. 1 are lines overlaying the data points, 
with fit parameters listed in Table I. Uncertainty bars are 
derived from counting statistics and represent one standard 
deviation. 



with the results of the fits to the data to Eq. 2 are shown 
in Fig. [3] and on a log-log scale in Fig. [4] the best fit 
parameters for Eq.'s 1-4 are shown in Table I. In fact, 
subtle differences between model fits are impossible to 
discern on the scale of Fig. [3] and are obfuscated due to 
different Tn values on plots like Fig. |4j but in the fol- 
lowing we elucidate the nuances of each fit. Fits were 
performed for reduced temperature within 0.1 of tN to 
0.01 of tN- No change in the width of the scattering as a 
function of temperature was observed. For the x = 0.022 
sample, an abnormal behavior was seen above the pre- 
ponderant transition temperature that was shown to be 
elastic by measuring with different energy resolutions of 
1 meV and 0.2 meV that showed identical behavior. This 
additional scattering was a consequence of the high sen- 
sitivity of the transition to cobalt doping; similar behav- 
ior was seen in a neutron diffraction study of potassium 
doped BaFe 2 As 2 P 
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5. (Color online) (color online) \ 2 surfaces from fitting 
y w j . Using the color map shown at the 



bottom of this figure, \ surfaces on a log scale illustrate 
the uniqueness of the extracted parameters as well as giving 
an idea of parameter correlation. In each subplot, a solid 
(blue) horizontal line illustrates the best fit ft value while a 
dashed (green) horizontal line shows the value for a relevant 
universality class (/3 = 0.326 of a 3D Ising model for x = 0.039, 
and P = 0.25 of a mean-field tricritical model for x = 0.022 
and x = 0.021.) 



Before delving into the physical relevance of the ex- 
tracted parameters, it is edifying to briefly examine the 
fits themselves. The quoted uncertainties are square 
roots of variances from the least squares algorithm. To 
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understand how the parameters interact during the fit 
procedure, the correlation coefficients are a useful, but 
not necessarily definitive metric,^ that suggest a strong 
connection between the proportionality constant, A, and 
the critical exponent, /3, and less so between other pa- 
rameters. A clear illustration of possible correlations as 
well as the goodness of fit are the x 2 surfaces near the so- 
lution, which we plot for f3 vs. Tjv and /3 vs. A for fits to 
Eq. 1 in Fig. [5j For a given data set, the four equations 
have similar x^maps, with log corrections (not shown) in 
eq.'s 3-4 systematically causing a shift to higher j3 values 
while retaining the basic shape of the minimum. 

IV. DISCUSSION 

In the present study, we have measured and analyzed 
the critical behaviors of the order parameter of cobalt 
doped barium '122' crystals with a precision approaching 
the experimental limit dictated by the inherent chemical 
inhomogeneities characteristic of such doped materials. 
We find that the shape of the onset of the magnetic tran- 
sition to be a function of the amount of doping, in a 
slightly more complicated way than previously hypothe- 
sized for these systems that suggested two distinct crit- 
ical exponents for '1111' and '122' materials depending 
upon whether and Tg were coincident or separated!^ 
This general trend of softening the transition with dop- 
ing is also qualitatively present in a study performed 
over a large range of doping values.^ Quantitatively, for 
the parent BaFe2As2, previous neutron diffraction ex- 
periments found /3 = 0.103±0.018, which is less than 
but near the 2D Ising value of j3 = 0.125P^ One pos- 
sibility for such a reduction in the effective exponent is 
the weakly first order nature of the magnetic transitions 
for doping values smaller than at the tricritical doping. 
Doped Ba(Feo.953Coo.o47)2As2 (x — 0.047) samples that 
are within the superconducting range were accurately 
modeled with (3 — 0.3, [6] which is near the expected value 
for a 3D Ising modelpH and we observe a virtually iden- 
tical value for x — 0.039, Table I. It is also worth noting 
that a similar j3- value was seen in nickel doped supercon- 
ducting samplesP^ 

In the iron pnictides, the tricritical point is of impor- 
tance because of its potential role in the onset of su- 
perconductivity. However, there is also a general inter- 
est in tricritical points from a fundamental viewpoint, 
stemming from the inception of the field in the study 
of 3 He- 4 He mixtures^ and quickly branching out to 
magnetic systems.^ Therefore, when additional systems 
come along that possess a tricritical point, it is exciting 
to test the validity of the theoretical models. In cobalt- 
doped BaFe2As2, we have studied the tricritical point in 
the (x,T) plane, which is expect ed to behave in a clas- 
sical manner such that f3 = 0.25.^^1 As we previously 
noted, logarithmic corrections to mean-field tricritical ex- 
ponents have been predicted by theorypS and such cor- 
rections have been applied to similar systems Nearly 



identical behavior is seen in the x = 0.021 and x = 0.022 
samples that are in the vicinity of the tricritical point. 
The best fit exponent for simple power law fits including 
the evident spread in Tjv due to chemical inhomogeneity 
(Eq. 2) is 0.21±0.01; this is somewhat less than the mean 
field tricritical value of 0.25, presumably due to the fact 
that our data do not probe the true asymptotic critical 
region. It is puzzling that the inclusion of a possible log- 
arithmic correction causes the values of beta extracted 
from the fits to increase dramatically, taking on physi- 
cally unrealistic values. For the x = 0.039 sample, where 
we expect a simple power law to describe the behavior of 
the order parameter well we find /3=0.30±0.01, close to, 
but somewhat less than the 3D Ising value of 0.326. This 
is consistent with the results of Wilson and coworkers^ 
who analyzed similar data in a large number of materials. 
It is possible that the reduced effective exponent reflects 
a residual effect of the nearby tricritical point. 

Finally, the moment values that we measure are in 
the expected range for the doping values measured,^ al- 
though the anomalously low value of the x = 0.022 mo- 
ment is unexpected. Previous systematic work has shown 
appreciable scatter in the moment values, and it is likely 
that strains in the sample contribute to this distribution. 
It is circumstantially evident from the additional features 
above TN that the x — 0.022 sample may have a larger 
internal strain. 



V. CONCLUSIONS 

In summary, we have measured the critical exponents 
of Ba(Fei_ a; Co a ;)2As2 around and above the tricritical 
point in the (x,T) plane, showing behavior consistent 
with mean- field predictions at the tricritical point. When 
above the tricritical point but below optimal doping, we 
find values of /3 « 0.30, consistent with previous results 
in a variety of materials but slightly below the expected 
asymptotic 3D Ising value of 0.326. As the presence of a 
tricritical point seems to be a common feature in many 
iron superconductor systems, it will be interesting to see 
if other systems show the same behavior. Finally, there 
should also be dramatic signatures of the tricritical point 
in both the heat capacity and the staggered susceptibility. 
These properties will be explored in future experiments. 
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